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Abstract 

Sampling and reconstruction of functions is a central tool in science. A key result is given by 
the sampling theorem for bandlimitcd functions attributed to Whittaker, Shannon, Nyquist, and 
Kotelnikov. We develop an analogous sampling theory for operators which we call bandlimited if 
their Kohn-Nirenberg symbols are bandlimited. We prove sampling theorems for such operators and 
show that they are extensions of the classical sampling theoremF] 



1 Introduction 

The classical sampling theorem for bandlimited functions states that a function whose Fourier transform 
is supported on an interval of length Q is completely characterized by samples taken at rate at least 
I/O per unit interval. That is, with J- denoting the Fourier transform 2 ! we have 



Theorem 1.1 For f G L 2 (R) with supp Tf C [-§, choose T with TQ < 1. Then 

\\{f(nT)}y {z) =T\\f\\ L i m , 
and f can be reconstructed by means of 

f( X ) = y f(nT) sin( :y (x - " )} 

neZ v 1 

with convergence in L 2 (M). 



Theorem 1.2 is an exemplary result from our sampling theory of operators. We choose a Hilbert- 
Schmidt operator H on L 2 (M) with kernel kh and Kohn-Nirenberg symbol an, that is o~h{x,D) = H 
in the sense of pseudodifferential operators |Hor791 Tay81| . Recall that Hilbert-Schmidt operators on 
L 2 (M) are exactly those bounded operators H with an £ L 2 (M 2 ) and corresponding norm of H. Let J- s 
denote the so-called symplectic Fourier transform 2 on L 2 (IR 2d ). 

Theorem 1.2 For H : L 2 (K) — > L 2 (R) Hilbert-Schmidt with supp T s a H C [0, T) x [-§,§) and T9, < 
1, we have 



<5fcT||i 2 (K) = T\\H\\ HS , 



1 201Q Mathematics Subject Classification. Primary 42B35, 94A20; Secondary 35S05, 47B35, 94A20. 
2 See Section [2] for basic notation used throughout this paper. 
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and H can be reconstructed by means of 



v^/ TT v^^ \/ r^-, s'm(irT(x — n)) 



with convergence in HS(L (. 



nez kez 



As shown in Section [3j Theorem 1 1 . 1 1 can be deduced from the general form of Theorem 1.2 which is 
stated below as Theorem 13,41 



The appearance of the sampling rate T in the description of the bandlimitation of the operator's 
Kohn-Nirenberg symbol reflects a fundamental difference between the sampling of operators and the 
sampling of functions. This phenomenon is illuminated in terms of operator identification by Theorem 



3.6 in [KP06 and Theorem 1.1 in [PW06a], results which we extend here by Theorems 4.6 and 4.7 
below. In fact, in the classical sampling theory, the bandlimitation of a function to a large interval 
can be compensated by choosing a correspondingly high sampling rate. In the here developed sampling 
theory for operators, only bandlimitations to sets of area less than or equal to one permit sampling and 
reconstruction. The bandlimitation to, for example, a rectangle of area 2 cannot be compensated by 
increasing the sampling rate, and, in fact, operators characterized by such a bandlimitation cannot be 
determined in a stable manner by the application of the operator to a single function or distribution, 
whether it is supported on a discrete set (which we shall refer to as sampling set below), or not. 



Theorems |4.6| and 4.7 in simple terms is Theorem 1.3 below. It can also be deduced from earlier 



operator identification results in [Pfa08a, PW06bj. As it is customary to define Paley- Wiener spaces 

PW(M) = {/ G L 2 (R d ) : supp Ff C M} 

to describe spaces of functions bandlimited to M C we introduce in this paper operator Paley- 
Wiener spaces 

OPW(M) = {H G HS(L 2 (R d )) : supp.FV.ff C M} 

to describe operators bandlimited to M C M. 2d . In short, PW(M) and OPW(M) are linked via the 
Kohn-Nirenberg correspondence [Fol89, KN65 . 

Theorem 1.3 Let (i(M) denote the Lebesgue measure of the set M C M 2 . 

1. For M compact with fJ>(M) < 1 exists T > 0, a bounded sequence {ck}, and A, B > with 

A\\H\\ HS < \\HY,c k S kT \\mR) < b\\h\\ hs , H g OPW(M). 

2. Let M be open with fj,(M) > 1, then, for all g G S'(R) and e > 0, exists H G OPW(M) with 

\\Hg\\ L 2 (m < e\\H\\ HS . 



The sampling theory developed here has roots in the work of Kozek, Pfander [KP06] and Pfander, 
Walnut [PW06aJ which addressed the identifiability of slowly time-varying operators, that is, of so- 
called underspread operators. Measurability or identifiability of a given operator class describes the 
property that all operators of that class can be distinguished by their action on a well chosen single 
function or distribution. The importance of operator identification and, therefore, operator sampling in 
engineering and science is illustrated by the following two examples. In case of information transmission, 
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complete knowledge of the communications channel operator at hand allows the transmitter to optimize 
its transmission strategy in order to transmit information close to channel capacity (see, for example, 
|Gol05| and references therein). In radar, simply speaking, a signal is send out and the goal is to 
determine the nature of reflecting objects from the received echo, that is, from the response to the radar 
channels input signal |Sko80| lKP06j . 

The best known operator identification example states that time-invariant operators are fully char- 
acterized by their response to a Dirac impulse. Kailath |Kai62j and later Bello |Bel69j investigated 
the identifiability of slowly time varying channels (operators) which are defined by the support size of 
their spreading functions, namely of the symplectic Fourier transform of the operators' Kohn-Nirenberg 
symbols. In both papers, conjectures were made that were then proven in [KP06], respectively [PW06aJ. 
The operator sampling Theorems 4.6 and 4.7 extend the main results in [KP061 IPW06a| . Shannon's 
sampling theorem, as well as the fact that time-invariant operators are identifiable by their impulse 
response (see Figure [I]) . 

The paper is structured as follows. Section [2] provides background on time-frequency analysis 
of functions and distributions, in particular on modulation spaces (Section 2.1), as well as on time- 
frequency analysis of pseudodifferential operators (Section |2.2[ ). In Section 2.3 we discuss boundedness 
of pseudodifferential operators on modulation spaces. In Sections [3] and [4] we state and prove our main 
results. Section [5] contains references to recent progress and open questions in the sampling theory for 
operators . 



2 Background 

L 2 (R d ) denotes the Hilbert space of complex valued, Lebesgue measurable functions on Euclidean space 
R d |Fol99j . The Fourier transformation T , respectively the symplectic Fourier transformation J- s , is 
the unitary operator T : L 2 (R d ) — > L 2 (R d ), f H> / = Tf, densely defined by 

7(7) = / f(x) e- 2 ™T* dx, f G L 1 (R d )nL 2 (R d ) , 

JR d 

respectively J 78 : L 2 (R 2d ) — > L 2 (R 2d ) with 

F s F(t,v)=ff F{x,i)e- 2 ^ u) ^ )] dxdi= ft F(x, £) e -^"-^) dx d£, FeL 1 (R 2d )nL 2 (R 2d ), 

where [•, •] denotes the symplectic form on R 2d . Throughout the paper, integration is with respect to 
the Lebesgue measure which we denote by fj,. 

The Fourier transform defines isomorphisms on the Frechet space of Schwartz functions 5(M rf ) and 
on its dual S ; (R d ) of tempered distributions (equipped with the weak-* topology). Note that 5'(IR rf ) 
contains constant functions, Dime's delta S : f h-> /(0), and weighted Shah distributions J2nez d c nO~kT, 
T G (M + ) d , with {c n } having at most polynomial growth. 

Similarly to the Fourier transformation, the time shift operatorTt, t G R d , given by Ttf{x) = f(x — t) 
and the modulation operator M w , w £ R d , M w f(x) = e 2mw ' x f(x), act as unitary operators on L 2 (R) 
and they are isomorphisms on 5(M rf ) and ^'(IR^). Note that M w is also called frequency shift operator 
since M w f = T w f. Further, we refer to vr(A) = 7r(t, v) = M v Tt for A = (t,v) E R 2d as time-frequency 
shift operator. Note that we have T o ir(t, v) = e 2mtv ir{y, —t) o T, that is, Tix{t, u)f = e 2mtu -n(v, —t)f 
for / G S'{R d ). 

The goal of operator identification is to select, for given spaces X and Y of functions or distributions 
defined on R d and a given space of linear operators Ti mapping X to Y, an element g G X which induces 
a continuous, open, and injective map <J> ff : H — > Y(R d ), H h-> Hg (see Figure [2]). 
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Figure 1: In the one dimensional case, the herein developed sampling theory for operators applies to any 
pseudodifferential operators whose Kohn-Nirenberg symbol is bandlimited to a compact set of Lebesgue 
measure less than one (for example, the blue region above). The results extend the classical sampling 



theorem described in Theorem 1.1 which is equivalent to the identifiability of operators whose Kohn- 
Nirenberg symbol is bandlimited to a segment of the frequency shift axis (red). Also, the fact that 
time-invariant operators with compactly supported impulse response can be identified from their action 
on the Dirac impulse is a special case of our results since the Kohn-Nirenberg symbols of time-invariant 
operators are bandlimited to the time shift axis (green). 
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Figure 2: Illustration of the operator identification and sampling problem. We seek an element g G X 
in the domain of the operator class Z which induces a map from Z into the range space Y which is 
continuous, open, and injective. If we can choose g = ^ ■ Cj5 Xj , then Z permits operator sampling. 

Definition 2.1 Let X be a set, Y a topological vector space, and Z a topological vector space of operators 
mapping X to Y . The space Z is identifiable by g G X if<& g :Z — > Y, H i— > Hg is continuous, open 
and injective. In the case that Y and T~L are normed spaces, this reads: there exist A,B>0 with 



A\\H\\ z <\\Hg\\ Y <B\\H\\ z , H £ Z . 



(1) 



If we can choose g E X = X(M. d ) of the form g = Cj5 Xj , Xj E M. d and Cj G C for j G Ij d , as identifier, 
then we say that Z mapping X to Y permits operator sampling and we call {xj} a set of sampling for 
Z with respective sampling weights {cj}. We refer to g as a sampling function for the operator class Z. 

In the following, we shall abbreviate norm equivalences as the one given in using the symbol >c. 
For example, ([!]) becomes 

\\H\\ z ^\\Hg\\ Y , H £ Z . 



We shall describe in Section 2.1 the distribution spaces and in Section 2.2 the pseudodifferential 



operator spaces considered here. Section [273] discusses boundedness of the considered pseudodifferential 
operators on modulation spaces. 



2.1 Modulation spaces 

To describe the full scope of operator sampling, we need to employ recent results in time-frequency 
analysis, in particular, we have to enter the realm of so-called modulation spaces. As Theorems 1.2 and 



1.3 indicate, all results presented include the special case of Hilbert-Schmidt operators and the Hilbert 



space of square integrable functions as range space, and we advise readers without significant expertise 
in time-frequency analysis to focus on this case during a first reading. 

Feichtinger introduced modulation spaces in |Fei81j . Modulation space theory was further developed 
by Feichtinger and Grochenig as special case of their coorbit theory [FG89J: for p being a square 



6 



G.E. PFANDER 



integrable unitary and irreducible representation of a locally compact group G on a Hilbert space H 
and Y being a Banach space of functions on G, we consider, for appropriate ip £ H, the so-called voice 
transform V v : H — > Y given by V ip f{x) = (f,p{x)<p), x £ G. Given an appropriate Banach space 
Gelfand triple X C H C X' , the coorbit space My consists of those f £ X' with ||/||my = H^/Hy < 00 
|FK98] . 

The special case of modulation spaces is based on the Schrodinger representation of the reduced 
Weyl-Heisenberg group. The corresponding voice transform simplifies to the short-time Fourier trans- 
form, that is, for any Schwartz class function consider 

V v f(X) = (f,7r(X)<p) = F{fT t Jp){v), A = (t, v) £ R 2d , 

which is well defined for any / £ S'(R d ) |Gro01] . Note that throughout this paper, dual pairings (•, •) 
are linear in the first component and antilinear in the second. Moreover, any choice of ip £ S{M. d ) \ {0} 
can be used to define modulations spaces (with equivalent norms), but as is customary, we shall choose 
a normalized Gaussian, namely (p{x) = g(x) = 2* e ^U^Ua ; x £ M. d . 

The role of the Banach space Y in coorbit space theory is attained in modulation space theory 
by weighted mixed L p spaces which we shall describe now. For a measurable function / on R d and 
p = {pi,... ,Pd), 1 < pi, ■ ■ . ,Pd < 00, we define mixed L p {R d ) spaces by finiteness of 

\\f\\Lv={j (•••(/ (j \f(xi,...,X d )\ Pl dxl) dx 2 ) ...dx d -i) dx d ) , 

with the usual adjustments if some p^ = 00 [BP61]. The mixed P(Z d ) spaces are defined accordingly. 

Note the sensitivity to the order of exponentiation and integration. For example, for f{x,y) = 1 if 
\x — y\ < 1 and f{x, y) = else, we have sup^ J \ f{x, y)\dy = 2 but J sup x \f{x, y)\dy = 00, that is, 
/ £ L 1,00 (1R 2 ) but g i L°°> l (R 2 ) where g{x,y) = f{y,x). 

A locally integrable function v : M. d — > M.q with 

v{x + y) < v(x) v{y), x,y £M. d , 

is called sub-multiplicative weight. For example, w s {x) = (1 + ||x||) s , s > 0, is a submultiplicative weight 
on M. d . If v is a submultiplicative weight and the locally integrable function w : R rf — > M^J" satisfies for 
some C > 

w{x + y) < C w{x) v{y), x, y £ M d , 

then w is a v -moderate weight function. The class of v -moderate weight functions on M d is denoted 
by M. v {M. d ). Note that for s < 0, for example, l®w s (x,!;) = (1 + ||£||) s is not submultiplicative, but 
l®w s is 1®W- S -moderate. If w is a v-moderate weight function with respect to some submultiplicative 
weight, then we simply say w is moderate. Note that for any moderate weight function on ~R d exists 
7, C > with ^e~ 7 '^"°° < w(x) < Ce 1 ^^ 00 (see Lemma 4.2 in |Gro07j ). A moderate weight function 
w on M rf is a subexponential weight function if there exists 7, C > and < /3 < 1 with 

£ e -7lWl£> < w ( x ) < Ce^ x ^. 

Weight functions on discrete groups such as Z d are defined accordingly. See |Gro07| for a thorough 
discussion on the role of weight functions in time-frequency analysis. 

Given a u-moderate weight function w, then the Banach space LS,(M rf ) is defined through finiteness 
of the norm H/H^p = The space Lw(M. d ) is shift invariant, shift operators are bounded on 

Lw(^ d ) but not isometric if w is not constant. Replacing M. with Z d , or with a full rank lattice 
A = A7j d , A £ ~R dxd invertible, either equipped with the counting measure gives the definition of 
/2,(Z d ), respectively Z£(A). If to is a moderate weight on then its restriction to A, which we denote 
by w, is moderate as well. 
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Definition 2.2 For p = (pi, . . . ,p d ) and q = (q\, . . . , q d ), 1 < Pk, Qk < oo, and w moderate on R 2d , we 
define modulation spaces by 

M™(R d ) = {/ G S'(R d ) : V 3 f G L£«(R M )} (2) 

\Fei81\ {GroOiy . The modulation space Mw q {R d ) is a shift invariant Banach space with norm = 
\\wV Q f\\ L v, q . Ifw = l, then we write MP>i(R d ) = M p ' q (R d ). If pi = ...= p d and (& = ... = q d then we 
abbreviate M£ 1,9l (R d ) = M i Pl '-' Pd) ' {qu -' qd) (R d ). 

Below we shall use the fact that replacing g with any other tp G 5(R rf ) \ {0} in ^ defines the identical 
space with an equivalent norm |Gro01| . Note that if p\ < P2, qi < <?2, and itfi > CW2 for some c > 0, 
then M p \' Ql embeds continuously in M P2 2 ,q2 , and, consequently, if W\ x then M£' 2 9 (IR d ) = M£f(R d ) 
with equivalent norms. 

The space M 1 ^ 1 (R d ) is the Feichtinger algebra, often also denoted by S (R d ), and M co ' ao (R d ) is its 
dual S' (R d ). In fact, in general we have M£' 9 (R d )' = M p ' f ^ (R d ) for 1 < p, q < oo and i + ^ = 1, 
- + ^■ = 1. Note that M 1( ^ )U)s (]R ) is also known as Bessel potential spaces, in particular L (R ) = 
M 2 ' 2 (R d ). 

To illustrate the chosen order of exponentiation and integration in the definition of the modulation 
space M p,q (R d ) for d > 1 and p / q, we state exemplary that / € M^'^ s ' ( ' 4 ' 5 ' ) (M a! ) if and only if 

/(/(/ ( / | (! + \Ai + 4)' V 9 f(ti, t 2 , u u u 2 ) | 2 dti) 3 dt 2 ) 5 di*) * di/ 2 < oo. 
Clearly f®g G M^^' (91,5a) (K M ) if and only if / G iW2i' 9l (K d ) and p G M£ 2 2 ' 92 (IR d ). In this case, 

\\f®9\\ M {vi,P2)Mi.n) = ii/iimsv^ii^Ims?,'* 2 - 

For compactly supported and bandlimited functions, modulation spaces reduce to weighted mixed 
L p (R d ) spaces. The following is a straight forward generalization of results in |Fei89| IOko09| . 

Lemma 2.3 Let p = (pi, . . . ,p d ) and q = (q\, . . . , q d ) with 1 < pk, qk < oo, let w = w\ ® W2 be a 
moderate weight function on R 2d , and suppose M C R d compact. Then 

1- ll/IU™ ~ \\f\\Ll 2 , f € S'(R d ), supp/ C M; 

2- II/IIm^^II/IIl^, f^S'(R d ), supp/CM. 

Modulation spaces allow for descriptions based on growth conditions of so-called Gabor coefficients 
[GroOl]. These descriptions rely on the following terminology. 

Definition 2.4 Let X be a Banach space, 1 < p±, . . . ,p d < oo, and let w be moderate on the full rank 
lattice A. 

1. {(?a}a6A Q X' is called l^-frame for X, if the analysis operator C^ gx y : X — > i^,(A), / h-> 
{(/, 5a)}asA > is well defined and 



x 



l|{(/,<7A>}lk, f€X. (3) 



2. {^aIagA Q X is called lw-Riesz basis inX, if the synthesis operator D^ gx y : C(A) — > X, {ca}aga ^ 
^2\C\g\, is well defined and 



E c a5a||x, {c A }G?P(A). (4) 
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In the classical Hilbert space setting X = X' = H and Z£(Z 2d ) = Z 2 (Z M ), the above entails the 
definition of Hilbert space frames and Riesz basis sequences. In the Hilbert space theory, condition Q 
implies that C{ gx y has a bounded left inverse, but in the general Banach space setting, ([3j alone does 
not guarantee the existence of such a left inverse. Therefore, the condition of a bounded left inverse Cj? 
is frequently included in the definition of frames for Banach spaces |Chr031 IGro914 IFZ98j . 

Note that for any 1 < p < oo, w moderate, /^,-Riesz bases form unconditional bases for their closed 
linear span. This follows directly from Q and Definition 12.3.1 and Lemma 12.3.6 in |Gro01| . 

For g G and A being a full rank lattice in R 2d , we set (g, A) = {Tr(\)g}\<=A- For w moderate 



on R 2d , set w\ = w(X). Results as Theorem |2.5| are important tools in modulation space theory, see, 



for example, Theorem 20 in [Gr604] or Theorem 6.11 in [Gr607j. 

Theorem 2.5 Let 1 < p, q < oo and let w be moderate on M. 2d . Let A be a full rank lattice in M. 2d and 
g G S(R d ). 

1. If(g,A) is a frame for L 2 (M. d ), then (g,A) is an l^ q -frame for M% q (R d ) . 

2. If (g,A) is a Riesz basis in L 2 (M. d ), then (<?, A) is an t p ~ q -Riesz basis in Mw q (R d ). 

Proof. 1. Let g G 5(IR d ) with (<?, A) is a frame for L 2 (M. d ). Let g generate the canonical dual frame of 
(g, A) of (g, A) [GroOlj . We have g G S{R d ) |Jan95j and conclude that both, C (gA) : M% g (R d ) l^ q (A) 
and £>(g,A) : l q, 9 W — > M% q (M. d ) are bounded operators. As Dq^oC^ is the identity on L 2 (R d ), 
we can use a density argument to obtain f^AjoC^^) is the identity on Mw 9 (R d )- We conclude that 
C(g,A) is bounded below. 

The proof of 2. follows similarly. □ 
2.2 Time— frequency analysis of pseudodifferential operators 

The framework of Hilbert-Schmidt operators suffices to develop the basics of our sampling theory for 
operators. But important operators such as convolution operators, multiplication operators, and even 
the identity are not compact and thereby fall outside the realm of Hilbert-Schmidt operators. Rather 
than focusing on operators with kernel in L 2 (M? d ), we shall consider kernels and symbols in modulation 
spaces. 

To formulate a widely applicable sampling theory for operators, we use the general correspondence 
of operators to distributional kernels given by the Schwartz kernel theorem (see, for example, [H6r07j). 

Theorem 2.6 For any linear and continuous operator H : S(M. d ) — > S'(M. d ) exists a unique kh G 
S / (M 2d ) with (Hf,g) = {k h J® 9 ), f,geS(R d ). 

Alternatively to kh, we can consider the so-called time-varying impulse response hu G S"(IR M ) of 
H : S{R d ) — > S'(R d ) which is formally given by 



h H (x, t) = k h (x, x - t), Hf(x) = j h H (x, t)f(x - t) 



dt. 



The Kohn-Nirenberg symbol an of an operator H : S(R d ) — > S ; (M. d ) is densely defined by an 
Ft-t^ha, that is, 



*h(x, = j kh(x, x-t) e~ 2 ^ dt, Hf{x) = f a H (x, /(£) , 



,2-irixi; 
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|Fol89| lKN65j . Note that the n-th order linear differential operator D : f t-t J2n=o an ( x )f^( x ) nas 
Kohn-Nirenberg symbol ar)(x,£) = Y2n=o a n{x)(2i:i^) n which is polynomial in £. Pseudodifferential op- 
erator classes, for example, those considered by Hormander, have symbols crjj which are not necessarily 
polynomial in £ but which satisfy corresponding polynomial growth conditions [H6r07 . 

Additionally, in time-frequency analysis and in communications engineering, the spreading function 
rjH is commonly used to describe H: 

m =T s a H , Hf(x) = Jj \ H (t,v)M v T t f(x)dtdu. (5) 

Equation §5§ can be validated weakly by first integrating with respect to x in 



(Hf, <p) = J JJ w(t, f)7r(t, v)f[x)<p{x) dtdv dx = ( VH , Vftp), f,^eS(M. d ), 

where Vfip(t,v) = (ip,n(t,u)f) is the short time Fourier transform defined above. Equation ^ illus- 
trates that support restrictions on r\n reflect limitations on the maximal time and frequency shifts which 
the input signals undergo: Hf is a continuous superposition of time-frequency shifted versions of / 
with weight function tjh [GP081 IKP061 lPW06aj . Moreover, as h(x, t) = J rj(t, p)e 2mvx dv, the condition 
supp77#(i, •) C [— |, |], t E R, excludes high frequencies and therefore rapid change of the time-varying 
impulse response h(x, t) as a function of x. In the time-invariant case, k(x, x — t) = h(x, t) = h{t) is, 
in fact, independent of x. These observations illuminate the role of support constraints on spreading 
functions in the analysis of slowly time-varying communications channels [Bel64, Zad52 . Additional 
aspects on the use of pseudodifferential operator calculus in communications can be found in [Str06] . 



2.3 Boundedness of pseudodifferential operators on modulation spaces 



Theorem 3.3 in Sectio n [3 provides the upper bound in ([!]) for Theorems 1.3 3.4, 3.6, 4.2 



follows from Theorem 2.7 



and 



4.6 



It 



which generalizes Theorem 4.2 in [Tof07] as well as results in |CG03| ICza03| 



GH99, GH04, Tac9U ITof04j where, generally, the case p% = q^ and p4 = q^ in the notation below was 
considered. Recall that p' denotes the conjugate exponent of 1 < p < oo, that is | + h = 1. 

Theorem 2.7 Assume 1 < p\, P2,P3,P4, Qi, Q2-, <?3, <74 < oo with p^ < qz,qA, 

1111,1111 

H < 1 1 and H < 1 1 . (6) 

P2 Pi P3 Pa (12 qi 93 <?4 

Let the moderate weight functions w,wi,W2 satisfy 

w{x,£,v,t) > c — -r- 

wi(t - x, £) 

with c > 0. Then, for some C > ; 

\\L a f\\ Mir * < C |k|| M , P3 , 93 ),(,4, P4 ) II/IUSV 91 ' / G M^(R d ), a € M^^(R 2d ), (7) 

consequently, L a : MS,\ m (R d ) — ► MS,l' q2 (R d ) is bounded for a £ M^ m) ' {qi ' Pi) (R 2d ) and L a denoting 
the operator corresponding to the Kohn-Nirenberg symbol a. 



Theorem 2.7 is a consequence of the following lemma. 
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Lemma 2.8 Assume 1 < Pi,P2,P3,P4, 9l, 92, 93, 94 < oo with p 3 < pi,p 2 ,P4, 93 < 91,92,94, 

1111 ,1111 

1 = 1 and 1 = 1 . 

Pi P2 P3 Pa 9i 92 93 94 

Let the moderate weight functions w,w±,W2 satisfy 

w(x, t, v, £) <w\(t- x, £,)w 2 (x, v + Q ■ 
Then for &(x, £) = g(x)g(x — t), we have 



1^/8)50(1-1) (x,t,v,£)w(x,t,v,£)\ P3 dx) P3 dt) Vi dn q3 dv 



PA . 93 



< 



, fGM^(R d ),gGM^( 



where f®go{\ 1 ) (x, t) = f(x - t)g(x). 
Proof. For g, f € we compute 



g(x) f(x'-t')e 



g(x — x)g(x — x— (t — t)) dt dx 



g{x')e- 2 ™' y Q{x'-x) / f(s)e- 27Ti( - x '- s ^Q(s-(x-t))dsdx' 



g(x')e- 2nix '( u+ ^Q(x'-x)dx' / f(s) e - 2 ^Q(s-(x-t)) ds 



= V a g(x,v+OV f(x-t,0- 
Assume 1 < Pi,P2,P3,Pi, 9i, 92, 93, 94 < 00. For w = 1 and w\ = W2 = 1 we have 



Pi v 53 



94 J_ 
94 



(/(/(/ \Vvf®go(\^)(x,t,v,0\ Pa dx) P3 dt) p *d£) q3 dv) 



FA . 93 



< 





11 w, er lUn ii^o, * + er ikn j P3 P4 #j 93 ^; 

23 ^94 . 93 94 J_ 



P4 93 . 94 
P'i Pi Jt | 93 



< 



llVg/^H^ 




\\V s f(;0 P3 \\L^\\V s g(;u + Cr\\Ln) P3 dz^ ' " :; 

II^S 

12 93 1 



93 
L r 2 



93 1 



93 
L s 2 




\V s g P3S1 dx) sl m dt 



22 93 . 1 
1 p 3 jA s 293 



M p 3 r 1^93''2 ||<7||M J >3 3 1>93 S 2 



To apply Young's inequality to obtain (13), we assume p4 > p% and choose ri,si > 1 with 

1 1 P3 
- + - = 1 + —. 

n si p 4 



Similarly, to obtain (14), we use 94 > 93 and choose V2, S2 > 1 with 

1 1 93 
- + - = 1 + -. 

T2 S 2 94 
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To conclude our proof of the unweighted case, we set p\ = p^fi, qi = 93^2, P2 = P3S1, and q 2 = </3S2- 
As all factors must be greater than or equal to one, we require p±,P2 > P3 and qi,q2 > 93- Moreover, 
(15) and (16) need to be satisfied, this holds if and only if ([8]) holds. The case that for some k, = 00 
or qk = 00 follows from making the usual adjustments. 



The weighted case follows by simply replacing V$G with w VgG in equations (12) till (13), and then 
replacing V s f and V a g by w\ V g f and w 2 V s g. This is justified by □ 
Proof of TheoremWA Let f,g€ S{R d ) and H with a H G M ( - P[hq3 ^ ( - q4 ' P4 \R 2d ) . Then 



\(Hf,g)\ 



h(x,t)f{x - t)dt g{x)dx\ = \{h H J®9°{\ o))| 
\(*H, Ft^{f®9° (J "o))| 



< 



\°~h\ 



M; 



1 1 w 



(17) 



where we applied Holder's inequality for weighted mixed L p -spaces to obtain (flT} |Grd01| . To obt am 
0, it suffices to show F t ^{f®go ( \ ) G j^W^M^ for / G M^' 91 and g G Note that 

replacing q by any other test function in ^ leads to a norm equivalent to || • ||m£' 9 j an d we choose to 
show that for ^ = J 7 ^^©, (5(x,^) = g(j;)g(a; — t), we have that 



1 -1 ' 
1 , 



\ L p' 3 ,q' 3 ,q' 4 ,p' 4 - 



is bounded by the left hand side in §W§ for / G Af^' 91 and 5 G M?^ 2 . Note that as 

I ±VvJ r t ^tf<2>go ( i ~ ) I (*> f» *) = I h v ®f®9° ( i ~o ) I ( x > *> 0. f.^e 



(18) 



the boundedness follows from an adjustment of the order of exponentiation and integration in (10). 
Minkowski's integral inequality, namely, 



\w(x,-)F(x,-)\ p dx 



< 



LP 



\w(x,y) F(x,y)\ g dy) dx 



implies that \\w(x,y)f(x, y)\\ij>,<i < \\w(y,x)f(y,x)\\Lip ii p < q. Hence, if q\ < p\ and q' 3 < p' 4 , then we 
can move the t-integral in between the x-integral and the ^-integral and obtain that (18) is bounded by 



the left hand side of (10). 



We now prepare to apply Lemma 2.8 Observe that if we assume 

Pi,Pl,P2 > P3, <?4> ?i> 4i > <?3, Pi > <7 4 , 



(19) 



and 



1 



+ 



1 



1 1 



, 1 1 

and 1 — T 

qi ? 2 



Pi P 2 P3 Pi 
that is, PA,p'i,P2 < P3 and g 4 , q[, q 2 < 53 and p 4 < g 3 , g 4 , and 



1 1 

- + - 

% 1i 



1 

Pi 

the latter being 



1 1 1,1 
+ 1 = 1 hi and hi 

P2 P3 Pi q\ 



1 

q2 



1 1 1 1 and 1 1 

pi P2 P3 Pi qi q2 



Q3 



1 1 

- + — 
53 qi 



+ 1 



1 

qi 
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and 



1111 

1 + — = — + — + — 
P2 Pi P3 Pi 



and 1 + 



92 



1 1 1 

— + — + —. 

91 93 94 



Hence, we obtain if ( 19 ) is satisfied. Note that for p < p and g < g we have M%] q embeds continuously 
in Mw ,q (see, for example, Theorem 12.2.2 in |Gro01| ). Hence, ([7]) remains true if we decrease Pi,P3,P4 
and gi, g3, q^, and/or increase p2 and q2- We conclude that ([7]) holds if ^ and j>4 < qs, q\ are satisfied. 

□ 



Remark 2.9 Note that for Hilbert-Schmidt operators, we have 

ll-f^ll-ffs = II^hIIl 2 = II^hIIl 2 = IIchIIl 2 = II^IIl 2 ) (20) 

a fact which is helpful to obtain norm inequalities of the form Q. But when considering modulation 
space norms for operator symbols, the chain of equalities p0[ ) fails to hold. For example, we have 

\(h H ,Tv(x,t,v,€)g)\ = \{a H ,ir{x,£,v,t)Q)\ = \{VH,n(t,v,£,x)g)\, 

but due to the implicitly given order of exponentiation and integration, 

II^H||M(P1.P2).(91,92) ^ ll CJ /i'llA/(Pl'<'2)>(<?l,P2) ^ \\VH || Af(P2.9l).(«2.Pl) ! H '. S (R ) > e>'(IR d ). 

Consequently, when defining a modulation space type norm on sets of pseudodifferential operators, 
one can base it on either hn, &h, or r]H, each choice leading to different operator spaces and norms. 
Lemma |2 . 8| gives a hint that it may be advantageous to define operator modulation spaces 
OMP 1 >P 2 > qi > q2 (L 2 (R d ),L 2 (R d )) on L 2 (R d ) through finiteness of the norm 

//* p p P2 91 *?2 1 

(/(/(/ \V£a H (x,t,Z,v)w(x,t,Z,v)\ Pl dx)^ dt)™ dt) 1 * dv)™ , 

where the symplectic short-time Fourier transform V s with respect to the window function g G 5(IR 2d ) 
is given by 

V B s F(x, t, £, i/) = T S (F ■ %Jq) (t, „), F G S'(R 2d ). 

This choice of order of exponentiation and integration respects arranging the time variables ahead of 
the frequency variables, while listing first the absolute time variable x and then the time-shift variable t, 
respectively, we first list the absolute frequency variable £ and then the frequency- shift variable v. More 
importantly, with this choice, we have 

||#/||mp2,92 < C\\H\\ MP3^, q3 ,u H G OM p ^^{M p ^{R d ),M p ^{R d )) 

for all 1 <Pi,P2,P3,P4, gi,?2j 93, 94 < °o satisfying 

For simplicity of terminology, we avoid the use of operator modulation spaces and symplectic short- 
time Fourier transforms in the following. Lemma [2 . 3| implies that this does not lead to a loss of generality 
in case of the here considered operator Paley- Wiener spaces. 



3 Sampling and reconstruction in operator Paley— Wiener spaces 

We introduce operator Paley- Wiener spaces. 
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Definition 3.1 For 1 < p, q < oo and a moderate weight w on R 2d , operator Paley-Wiener spaces are 
given by 

OPW™{M) = {H : S(R d ) — > S'{R d ) : suppT s a H Q M and a H € L™{R 2d )}. 

OPWw ,q (M) is a Banach space with norm \\H\\oPWS,' q = II^IIlS,' 9 • If w = ^ and p = q = 2 then we 
simply write OPW(M) = {H G HS{L 2 (R d )) : supp T s a H C M] . 



Note that, as illustrated in Corollary |3.7| and Example |3.8| below, it is appropriate to choose 
OPWw°°(M), respectively OPW£? ,q (M), when considering multiplication respectively convolution op- 
erators. Moreover, observe that OPW£? ,CO (M) consists of all operators in the weighted Sjostrand class 
with Kohn-Nirenberg symbol bandlimited to M [Gro06, Sj694, Sj695, ,Str06j. 

Remark 3.2 In [H6r07], Hormander considers pseudodifferential operators with Kohn-Nirenberg sym- 
bol in 

S™ = [ a G C°°(R 2 d) : \d£dgtr(x,Z)\ < C a ^(l + ||^|| 2 )«-p(«i+-+««0+^i+-+^) j a, /3 G Nq} 

where m G R, < p < 1, and < 8 < 1. Clearly, if suppler C M and cr G S 1 ^, then L CT G 
O^^(M) if s < -m and L a G OP^(M) if (m + s)q < -1. 

Theorem 3.3 Zei 1 < p, q < oo and to moderate on R 2 . For M compact exists C > with 

\\Hf\\ M ™ <C\\<r H \\^ ||/||w-.ac, G OPW™(M), / G M 00 ' 00 ^). 

Consequently, any H G OPW^,' q (M) extends to a bounded operator mapping M 00 ' 00 (R d ) to Mw ,q (R d ). 

Proof. Set u(x,£,v,t) = w(x,£ + z^) and choose G with suppt/? C [-l,l] d . Then we use 

Lemma 



2.3 



and supp V^pcra Q ^ 2d x M + [—1, l] 2d , hence, uj x w®l on supp V^ipCXH, to obtain 



An application of Theorem 2.7 with f>i = gi = oo, that is p' x = q[ = 1, <pi = P3 = P> 12 = Q3 = Q, and 
P4, <?4 = 1 concludes the proof. □ 
In the following, we set Q T = [0,Ti)x . . . x[0,T d ) for T = (Ti,...,T d ) G and R n = 

[-% ^)x . . . x[-^, ^) for = (fix, . . . , n d ) e (R + ) d . 

Theorem 3.4 Let 1 < p, q < oo and Zet w; = it>i<8)it;2 oe moderate on R 2d . Let T, Q G satisfy 
TmQm < 1, to = 1, . . . , d. Let A = TiZ x . . . x and choose s G M 1,l (R d ) with supp s C Ri/t and 
se 1 on Rq. Then 

\\Hj2h\\M^-\\H\\oPW^, HeOPW™(Q T xR n ), (21) 
AeA 

and any H G OPWw' q {QT^Ro) can be reconstructed by means of 

k h {x +t, x ) = X Q T (t) Y,( H Y1 M(* + A ) s( - x - A )- ( 22 ) 

AeA A'eA 

with convergence in OPWw 9 (R 2d ) for 1 < p, q < oo and weak-* convergence else. 
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Proof. We shall show (22). The norm equivalence (21) can be shown by adopting the steps of the 



proof of Theorem 4.6 



For A = T{L x . . . x T^Z, we consider the Zak transform given by 



AeA 



Note (tf J] 5 A ,)(x) = (M^v), ^ <M = E 

kh(x,X') = y~] hff(x,x — X'). We consider first 



A'eA 



A'eA 



A'eA 



A'eA 



hu G M 1 ' 1 (M. 2d ) and use the Tonelli-Fubini Theorem and the Poisson Summation Formula |Gro01j . 
page 250, to obtain for (t, v) G Q T j_ 



2ni\v 



A'eA 



Y E h H (t- \,t- \- \')e 
AeA A'eA 

f ~ X ~ A '' v')e 2m(t ~ Xy du' e 2mX » 

AeA A'eA 

E E / Wit - A - A', l/' + I/ ) e »rt((*-A)(H-0+^) 
AeA A'eA 



i/')e w ' e " 2 ™ Al/ ' eft/' 



A"eA AeA 

^ E E »7H(t-A"^ + A)e 
A"eA AeA^ 



2mt\ 



A"eA AeA x 



where A 1 = {A G l d : e 2?riAA ' = 1 for all A' G A} = jrZ x . . . x is the dual lattice of A. 
This leads directly to (22) since 



/ XQ T (t) E(^ E^')(* + A ) s ^-A)e- 27ri - 
AeA A'eA 

- XQt(*)£ (^E 5a ')(* + A ) fs(x-X)e 



-2irivx 



dx 



AeA A'eA 
XQ T (i)£ (ff£»(t + A) e" 2 ^^) 



AeA A'eA 



.(t)S(i/) (Za(^EM) (*. 
V A'eA / 



= r/ H (t,i/)e^ 
We can apply Lemma 



h H {x,t) e- 2 ™( x -Vdx= h H (x + t,t) 



-2-kivx 



dx. 



2.3 



to show that ||/f|| OPW -p,9 X ||^if|| M (_ P ,i),(i,?), w(x,t, v,£) = w(x,£), and 



validity of ^ for h H G mJ^ ),(1 ' 9) (K m ) follows then from the density of M~' L (M? d ) in Mif ±MX,yj ( 
In case of p = oo or g = oo it follows from weak-* density. 



~w 



□ 



Remark 3.5 If p, q = 2, then we can alternatively choose s = XRn £ M 2,2 (M d ) = L 2 (M d ). This allows 
us to replace the inequality T m Vl m < 1 by T m Q m < 1 in the hypothesis of Theorem |3.4| 
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Note that Theorem 3.4 and its proof generalize trivially to the following setting. 

Theorem 3.6 Let 1 < p, q < oo and w = w\®W2 be moderate on ~R 2d . Let A,BQ~R d be bounded, and 
let A be a lattice such that A is contained in a fundamental domain of A and for some e > 0, B + [—e, e) d 
is contained in a bounded fundamental domain of A 1 - = {A E R rf : e 2mXX = 1 for all X' E A} . Choose 
s E M 1 ' 1 ^) with supps C B + [-e, e) d ands = l on B. Then 

ll# E^IU- ~ ll#llcw™, H E OPW^(AxB), 
AeA 

and any H E OPWw ,q (AxB) can be reconstructed by means of 

K H { X +t,t) = X A(t) Y,( H Y1 Sx ') (* + A ) s ( x - A )- 

AeA A'eA 

with convergence in OPWw' q (AxB) for 1 < p, q < oo and weak-* convergence else. 

Considering OPW P '°°([0, T)®[— ^, ^)), we obtain the classical sampling theorem as corollary to 
Theorem I3.4L 



Corollary 3.7 For m E L P (R), 1 < p < oo ; wif/i suppm C and T with TQ < 1 choose 

s E M 1,:L (IR) toitt supps C [—§,§) and seI on [-^, i). T/ien 

||m||x,p x || {m(kT)} \\i P (23) 

and 

m(x) = m(kT) s(x — kT). 

Proof. For m E with suppm C [— ^, we define the multiplication operator M formally by 

M : / — >m-f,fe S(R). We have 

M/(x) = m(a?)/(a?) = J m{x)5 (t)f(x-t)dt = J J 5 (t)m(u)e 2nixu f(x - t) dt dv . 
Hence, 6 ®m = r) M = T s o u , and, picking any T with TQ < 1, we conclude M E OPW p >°° ([0, T)x [-§, §)) 



Theorem 3.4 implies that -ff and therefore m is fully recoverable from M Ylk&z ^kT = Efcez m (kT)5kT, 



in fact, the reconstruction formula ([22j) reduces then to the classical reconstruction formula for functions: 
m(x)5 (t) = m(x + t)S (t) = n M {x + t,x) Th =^ X[0 T) (t) ^ (M y~] 8 kT ) (t + nT) s(x- nT) 

= X[o,T)(*) m ( kT )$kT{t + nT)s{x - nT) 



X[o,T)(t) E m ( kT ) S o(t - (kT - nT))s(x - nT) 

' 0, if t$ [0, T) or t £ 7LT; 

Enez Efeez m(A;T)5 ((n - fc)T)s(a; - nT), 

= J2 k £ Z m(kT)s(x - kT), if t = 0. 

5 {t)^2m(kT)s(x - kT) . 

fcez 
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The norm equivalence in ( 23 ) is obtained by verifying that 



|m|| LP x ||M|| iw.°° X ||M^J<5 nT || M p,cx. = || ^2m(nT)5 nT \\ M p,o° x \\{m(nT)} n \\ p , 

me L P (R), suppmC [-^,^]. 



□ 



In addition to the application of Theorem 3.4 to multiplication operators, we consider now 
OPW°°'P([0,T)®[-%, §)) for convolution operators. 

Example 3.8 Time invariant operators are convolution operators, that is, 

Hf(x) = h * f(x) = Jh(x- s)f(s) ds. 

Such operators represent the classical example of operator identification/sampling namely, as H5q(x) = 
h(x), H5q determines h and therefore H completely. In the framework of operator sampling, we 
consider h £ L P (R) with supp/i C [0,T]. We have r) H {t,v) = F s an{t,v) = h(t)5o(u) and H € 
OPW°°' P [[0, T) x [— |L, 4^))- Moreover, with appropriate s we may obtain H5q = h from (22), as 



h(t) = h H (x,t) = h H (x + t,t) Th ^ m X[0!T) (t)Y, {HY,5kT){t + nT)s(x-nT) 

= X[o,T) (t) (kT -nT))s(x-nT) = J2Yl *1°> T ) WW ~ ( kT ~ nT ^ s ( x ~ nT ) 

= H5 (t)s(x - nT) = H5 (t) ^ s{x - nT) = F5 (i) ^ s(|)e 27rirf = H8 (t)s(0) = H5 (t) . 

The distributional spreading support of a time invariant operator is also indicated in Figure [TJ 

4 Necessary and sufficient conditions for operator sampling and iden- 
tification 

The aim of this section is to show that the applicability of sampling methods for operators depends 
solemnly on the size of the spreading support set M, that is, on the Jordan content of M (see Defi- 



nition 4.4 below). Our main result in this section, namely Theorem 4.6, though, only covers the case 
d = 1, that is, operators H : S(R) — > S'(R). Possible means for generalizing Theorem 4.6 to operators 
H : S(R d ) — ► S'(R d ) are briefly discussed in Section || 

Before recalling the definition of Jordan domains and some of their properties, and before stating 



and proving Theorems 4.6 and |4.7[ we will use a geometric approach to obtain a sufficient condition 
for the identifiability of OPWl' p {M) if M = A(Q T xR n ) + (t ,i/ ) C R 2d , T, Q G (R+) d , T m Q m < 1, 
m = 1, . . . ,d, and A is a so-called symplectic matrix. Theorem |4.2| below generalizes Theorem 5.4 in 
[KP06j . 

Definition 4.1 The symplectic group Sp(d,R) consists of those matrices A £ SL(2d,R) = {A £ 
R 2dx2d . detyl = ^ with A * ( 0^ -I d \ A = ^0 -I d \ ^ where j d is the dxd identity matrix. 

Note that A £ Sp(d, R) if and only if [A(x, £) T , A(x', £')T] = [(x, £), (x 1 , £')] where [•, •] is the symplectic 
form defined in Section [2] 
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Theorem 4.2 Let A E Sp(d,R), (in^o) G H^ 2rf > 1 < P < oo, and Zei w be a moderate weight on R 2d 
with w(A(x,£) T ) < w(x,£). Then 

1. OPWS,' p (M) mapping M 0O ' co (R d ) to M p ' p (R d ) is identifiable if and only if OPWS,' p (AM+(t , v )) 
mapping M°°'°°(R d ) to M£j' p (R d ) is identifiable, and, consequently, 

2. for T, ft E (R + ) d with T m n m < 1, m = 1, . . . , d, we have OPW% p (A(Q T xR n ) + (to, uq)) mapping 
M°°>°°(R d ) to M p ' p (R d ) is identifiable. 



The proof of Theorem 4.2 is based on the representation theory of the Weyl-Heisenberg group. Here, we 
only outline the proof, the interested reader can import details from Section 5 in [KP06J or [Fol89"t lGro01| . 

Proof. We will obtain the identifiability OPWS,' p (AM) with A E Sp(2d, R) from the identifiability of 
OPWw P (M) by using the canonical correspondence between elements in OPW p ' p (AM) and elements in 
OPWw P (M) which is given by a coordinate transformation in the spreading domain R 2d D M, AM. In 
fact, Theorem 5.3 in |KP06| recalls that for A E Sp(d, R), there exists a unitary operators Oa on L 2 (R d ) 
with ir(A(t, u)) = OA^{t,v)OA*, t, v E R d . Such operators Oa, A E Sp(d,R) are called metaplectic 
operators, and they are intertwining operators for representations of the reduced Weyl-Heisenberg group 
that are unitarily equivalent to the Schrodinger representation [Fol89| iGroOlj . Metaplectic operators 
are finite compositions of the Fourier transform, multiplication operators with multiplier e~ nix with 
C selfadjoint, and normalized dilations / i— > |det-D|5 f(Dx), D invertible. They extend, respectively 
restrict, to isomorphisms on M p,p (R d ), 1 < p < oo, if w(A(x,£) T ) < w(x,£,) (see Theorem 7.4 in 
[FG92]). 

The following formal calculations of operator valued integrals can be justified weakly for all H E 
OPW p,p (AM). A similar computation can be made for H E OPW p ' p (M + (to,fo)), combining both 
leads Theorem |4.2[ We compute 



H = JJ r]H(t,v)M v T t dtdv = JJ rj H {t,u) <n(t,v)dtdv 

= J J n H (A(t, v) T ) Tr(A(t, v) T ) dtdv = J J n H (A(t, v) T ) A n{t, u)O a * dt dv 
= A [[ VH A (t, v ) v ) dt dv °a* = A H a O a \ 



where t] Ha = i] H oA and H A E OPW% p (M). The identifiability of OPW$ p (M) with identifier f M E 
M 00 ' 00 (IR d ) leads therefore to the identifiability of OPW% p {AM) with identifier f AM = A fu E 
M 00 ' 00 ^). In fact, we have 

||-H"/am||mS' p = W H °AfM\\MP' p ~ \\OA*HO A fM\\MZ P 

= \\HAfM\\ M ^-\WHj L ^^hH\\ LPwP = \\H\\ OPW ,, P , HGOPW p ' p (AM). 

□ 



Remark 4.3 Theorem 4.2 is not an operator sampling result per se as not necessarily all Oa map 



discretely supported distributions to discretely supported distributions. But Theorem |4.2| can be used 
to show that OPW p ' p (M) permits operator sampling by showing that 

1. M C AM + (to, ^o), 

2. A E Sp(d,R), 



3. w(A(x,0)<w(x,0, 
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4. OPWw P (M) permits operator sampling with sampling set {xj} and weights {cj}, and 

5. 0* A J2 c j$Xj is discretely supported. 



Note also that the restriction to p = q in Theorem 4.2 is necessary, as, for example, the Fourier transform 
is not an isomorphism on M p ' q whenever p ^ q. 



Theorem 4.2 relies on arguments based on symplectic geometry on phase space. As discussed above, 
Theorems 



4.6 



and 4.7 give a characterization for the identifiability of operators H : «S(K) — > <S'(M) 



which does not rely on any geometric properties. 
Definition 4.4 For K,LgN set R k ,l = [0, ^) x [0, f ) and 



J 

Uk,l = { |J (Rk,l + (|, ¥)) : k iiPi g J e n}. 



27ie inner content, respectively outer content, of a bounded set MCI 2 is 

vol - (M) = sup{//(J7) : UQM and U € /or some € N} , (24) 

respectively 

vol+(M) = inf{/i(C/) : COM and C7 G /or some K,L G N}. (25) 

Clearly, we have voP(M) < vol + (M). I/vol~(M) = vol + (M), i/ien we say i/iaf M is a Jordan 
domain with Jordan content vol(M) = vol~(M) = vol + (M). 

We collect some well known and useful facts on Jordan domains to illustrate their generality [Fol99j . 

Proposition 4.5 Let MCI 2 . 

1. If M is a Jordan domain, then M is Lebesgue measurable with /x(M) = vol(M). 

2. If M is Lebesgue measurable and bounded and its boundary dM is a Lebesgue zero set, that is, 
H(dM) = 0, then M is a Jordan domain. 

3. If M is open, then vol"~(M) = fi(M) and if M is compact, then vol + (M) = n(M). 

4- IfVQNis unbounded, then replacing the quantifier "for some LcN" with 'for some L £ V" in 



(24) and in (25) leads to equivalent definitions of inner and outer Jordan content. 



The second main result of this paper has been stated in simple form as Theorem 1.3 part 1, in 
Section [T] It also generalizes Theorem 1.1 in |PW06a| . 

Theorem 4.6 For 1 < p, q < 00 and w = w\ W2 moderate, the class OPW%} q {M) mapping M°° ,00 (M) 
to Mw' q (^) permits operator sampling z/vol + (M) < 1. In fact, i/vol + (M) < 1, then there exists L > 
and a periodic sequence {c n } such that 

\\Hj2^n\\ M P, q >,\\H\\ OPW ^, HGOPW^(M). (26) 



Sampling of Operators 



19 



Theorem 4.6 is complemented by Theorem 4.7 which generalizes Theorem 1.1 in |PW06aJ and Theorem 



5.2, part 2, in |PW06bj . 

Theorem 4.7 Let 1 < p, q < oo and w subexponential. The class OPWw q (M) mapping M 00 ' 00 (]R) to 
M£' 9 (M) is not identifiable if vol' (M) > 1, that is, for all f G M 00 ' 00 we have 

\\Hf\\ M ™ * \\H\\opw^ H e OPW^(M) . 



Theorem |4.6| is proven below. Subsequently, we outline the proof of Theorem 4.7 which employs 
elements of the proof of Theorem 1.1 in pPW06a] and Theorem 3.13 in |Pfa08b . 

4.1 Proof of Theorem I3~6l 

The following observations are special cases of Theorem |4.2[ They will be used in the following to reduce 
notational complexity. 

Proposition 4.8 Let 1 < p, q < oo and let w be a moderate weight on M?. 

14 



1. OPW% q {M) is identifiable by f if and only ifOPW% q (M- 1 / a °J ) is identifiable by D a f : 
f(ax). 

2. OPW£,' q (M) is identifiable by f if and only ifOPWS,' q (M + A) is identifiable by vr(A)/. 



Proof. We shall proof 1., the proof of 2. follows similarly. For H G OPW%J q {M- °J ), define 
H a G OPW% q {M) by rj Ha = Vnof 1 ^ ). Then a Ha = <J H o f V° °) as well. We compute formally 



{HD a f)(D = Ja H ^-,Oe 2 ^LXf(0^ = lJa H ^Oe 2m ^f(l)dC 
= J a H {% aO e 2 ™Z /(£) d£ = H a f(x). 
Using standard density arguments, we conclude that 

\\HDJ\\ M ^\\D,HD a f\\ k ^ H G OPW™{M- (V« °) ) 



□ 

Assume now that vol + (M) < 1. Applying Proposition 4.8, we assume, without loss of generality, 



that for 5 > sufficiently small, M + [-§, |) 2 C [0, 1) x M+. We choose K, L G N with L prime so that 
the following conditions are satisfied for some < e < 5 and M e = M + [— §, §) 2 

1. vol + (M e ) < 1, 

2. M e C [0,1) x [0,K), 

3. L > K, 

L-l 

4. M e C U M = (J (Rk,l + ¥)) , m^nj G Z, where i?^ L = [0, i) x [0, f ) and (m^nj) / 

i=o 



(mf,nj>) if j / j'. 
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Note that 1 = vol(U M )- 

The following result from [LPW05] is a key component of our proof of Theorem 4.6 In fact, if the 
restriction to L prime below could be weakened, then we would obtain a generalization of Theorem 4.6 
to higher dimensions (see Section [5]). 

Theorem 4.9 For c G C define 7i(k,£)c by (ir(k,£)c)j = Cj_ k e ^, k,£ = 0, ...,L - 1, where 
j — k is understood modulus L. If L is prime, then for almost every c G C L , the vectors in Q c = 
{Tr(k, £)c}k,e=o,...,L-i are ^ n general linear position, that is, any matrix composed of L vectors of Q c is 
invertible. 



Remark 4.10 Theorem 14.91 can be reformulated as a matrix identification result with identifier c 
[KPR08]. The use of algorithms based on basis pursuit to determine a matrix M from Mc efficiently is 
discussed in [PRT08L IPR09] . 



We now choose as c 6 Z°°(Z) the periodic extension of a vector (cq, . . . ,cl_i) which satisfies the 
conclusions of Theorem |4.9[ In the following, we shall show that kh can be recovered from Hg with 

9 = E»^el^(l). 

For simplicity, we shall assume first that kh G M 1 ' 1 (M 2 ). This additional assumption implies that 
for g G AI°° ,00 (M.), we have .Kg G M 1,1 (M.) [PW06bJ. This enables us to switch the order of integration in 
many of the following computations. After the necessary computations are completed, we shall extend 
our result using standard density arguments. 

Choose nonnegative (p G C%°(M) with J <p(x) dx = 1 and supp^J C [— |, |). We shall consider the 
case 1 < p < oo only, the case p = oo requires only the usual adjustments. 

Note that 



\V v Hg(t+%,v 



J2 c k h(x, x + %) e~ 2 ™ x v(x - (t + £ ) dx 
fcez 

J2c k h(x + %,x + § + ±) e~ 2mu % ir(t, uMx) dx 



E 

fcez 



c k h(x + %,x+ £ + ^)7r(i,z/)ip(x)cfo: 



Set u> = SmeZ^Cir) XrmL (m+i)L -. and observe that w x u>i. Then 



\\Hg\\ 



n + l 



|^#«7[ll#«= / " \V v Hg(t,u) Wl (t)\ P ) P dt 

n& K 



%)\ p Vdt 



V/ \V (p Hg(t + %,u)w 1 (t + 
Y, / '\V v IIg(t + %,vMt + %)\'y dt 
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We set ip = Xj_e K + i) * f an d observe that ip{u')T^ip{v') = T u ip(v') for oj G [0, ^), a fact that will 
be used to drop ip in (31) below. We compute for t G [0, ^) 

^Kt + f)| p |^(t + f^)f 

= ^H* + ^)I P | / E Cfc + + # + |)vr(t, u)ip(x) dx 

EEm^)I p | / E c ;-*M*+ 



, x + -^)ir(t, is)tp(x) dx 



j=0 m£Z 
L-l 



fcez 



Ell E / E c ^- fc/iH ( x+ ^^jr 2 '^ ^)^(t,i , )f(x)dx 



j=o mez " fcez 



(27) 



L-l 

E 

i=o 



fcez mez 



,mL+i 



ir(t, v)v{x) il>(i/) tt(x', g)<p{v') dv' dx" p 



1 Q w i 



(28) 



where we used that M mL+ 3 ijj is an /--Riesz basis in the Banach space Mj^ (R) to obtain (27), that 
is, we used 



E 



|| a m M_ m L± 1 4>\\ M P,P = || ^ a m V tp M_rnL± i ip\\ L P ii 

m£Z m£Z m£Z 



Note that 



■2niu 



I mL+j 
K 



m6Z 



E / 

m£Z 



E^+l^' + f )e 



(29) 



2 w i(x+^)(i/+^) e -2 W »(|,i/+^) e ari^ 



(30) 



te2 



We have applied the Poisson Summation Formula to obtain (29). Moreover, we chose r}nix! ,u' 
rj H {^ >')e Wi/ in d30b. 
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After substituting (30) into (28), we integrate with respect to t on [0, h) to obtain 



[ K J2 \Mt+%,")\ P KHg(t+%,v)\ p dt 



j=0 u feG2 



L-l „J_ 



7r(i, -0(1/') 7r(^',a;')</p(i/') dz/' <ix|| rP (it 



E jf // 1 / / 5 § c '-' ^""^ + + 

e" 2 "^ 1 '' 7r(t, i/)p(x) ^V^'M^)^' ™iCO T dx ' ^' dt 



i=o 



L-l „ J_, 



j=o Jo J Jo " fcez te2 



- w £ r^(x + 4, i/' + ) n(t, V(^) 7r(f, x')^(^) dv' dx wi(x')\ p d? dx' dt 



E^/flEE,-*^ 

e " w £ ^(s + + ir(t, u)ip{x) 7r{C,x')^(u') dx dv' w 1 (x')\ p d? dx' dt (31) 

././ 



E/ 7/ IEE<^ ! 



V^m(t + ±,Z' + f,v,x' + ±) e 2m ^ e 2 -f ^(x') | p d£' dx' dt 



j=0 f=0 



J Jo 



V^mit + i ^,v,x'+ k £) e 2 ^ v e 2mx>l£ ^ w x {x') \ p d? dx' dt (32) 

j?rrf E I w»(* + + + k 4) ^ v ^0*' + k 4)\ p d? dx' dt 

j'=0 

K 



0JJ0 Jl=0 



E/ // Iw^+^.e'+^^^o^i^or^'^'dt 

; . /=0 Jo JJO 

r (*:,/+!) (i ,-/+!) 

E /V \V^<prf(t,g,u,a/)wi(x')\ p dg dx' dt. 



j'=0" K 
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To obtain (32) we used that V v ^ v rj C [0, l)x[0,K). Moreover, we used 

,/A 



-2niu' 



K m(x + f , v' + Zf ) 7r(t, v)yj(x) vr(f , x')^') dx dz/ 



-2-kiv 



If ^(x + £ , v' + /£ ) e- 2 ™>(i - t) e" 2 ™ " <p(i/ - £') dx du' 
m(x, v') e- 2m ^ x -^ip(x - (t + I)) eT 27 ^'- 1 ? )F e -2*<*V-^) ^ _ g + IK^ dxdv , 



k 



Using the fact that replacing now <p (g) (p by any other test functions leads to equivalent norms of the 
modulation space at hand, we obtain for real valued g € S^M 2 ), g(t, u) = 1 for [—1, 2)x [-K, 2K), 



\\Hg\\ M ™ ~ \\VtpHg\lLP,, 



n+l 
K 



1 

\P\" 



L <V +1) 



> 



> 



J2 /v 

j'=0 J -K 



I \V v Hg(u,u)w 1 (u)\ 

ti6Z t 

Wpgnpfjiu, £, is, x)wi (x) \ p d^dxduV 




IQJO 
•lrK 



WO 
•l r K 



| Vtp^tp^u, £, ^, x)wi (x) | p d£ dx du J P 



| V^7/(n, £, i/, x)wi (x) \ p d^dxduy 



>oJo 

-lrK, 



I X[o,i) (u)x[o,i<) (£)V Q rj{u, £, u, x)wi (x) \ p d£,dxdu) P 



£<(!/) 




wo 



X[o,i) («)X[o,K) (O^viv, x )wi (x) \ p d£dx du 



r K \ 1 

K J \T s ?j(u,x)w 1 (x)\ p dx^ r 

(j\\ rP-9 X (T rP,<7 . 



z/)u>i(x)| p dx) P 



(33) 



To obtain (33), we apply a mixed L p -space version of Young's inequality for convolutions, namely, 



we use that for £>(x,£) = e 2mx ^g(x,^) € L 1 (M), we have cr(x,£) = g*a and a(x,£) = g*a (see Theorem 
11.1, [(iroOlp . 



4.2 Outline of the proof of Theorem 4.7 



We omit detailed computations as they would closely resemble computations carried out in [KP06, 
Pfa08bl IPW06bl IPW06a] . For the interested reader, we suggest to use (PW06aj as a companion when 
filling in detail. 

We shall show that for a measurable subset M with vol~(M) > 1, the operator class OPW p ' q (M) 
is not identifiable. In detail, we shall show that for every g £ M 00,00 (]R), the operator 

<5> g : OPW™(M) — ► M™(R), H h-> Hg , 

is not bounded below, that is, there exists no c > for which we have ||.ffg|| M p>9 > c||<7#|| i P,3 for all 
H G OPW p,q (M). 
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OPW£,' q (M) 



{^} 



C, 



(fl,«Z 2d ) 



D 



{Pj} 



C, 



(g,uZ 2d ) 



M 



M 



Figure 3: Sketch of the proof of Theorem 4.7 We choose a structured operator family {Pj} Q 
OPWw' q (M) so that the corresponding synthesis map D^py : {cj} — > E c j-^j nas a bounded left 
inverse. Note that C( 0j u ^2d) has a bounded left inverse for u < 1 as well. Theorem |4.13| shows that for 
any g G M°°(M), the composition M = C, 



that (j) g : OPW% q {M) 



7 2d^O 



)D{p.y has no bounded left inverses. This implies 



also has no bounded left inverses. 



To this end, choose K, L G N and Vm = \jf=o {^-K,L + > m j? n j e where -R/^l = 

[0, x [0, ^) and where (mj,rij) / (rrij',njt) if j ^ j', such that Vm ^ M and voI(Vm) > 1- It is 
sufficient to show that OPW% q (V M ) is not identifiable as OPW$ q {V M ) C OPVC(M). 

The proof of Theorem |4.7| is also sketched in Figure |3j The proof is based on extensions to results 
from [Pfa08b, PW06aJ which are stated below and which concern the construction of the operator family 
{Pj} in Figure [3} 

Lemma 4.11 Let P : S(R) -> S'(R). For p, r G E and to, £ G E ; de/me P = M w T p - r PT r M^ w : 
S(R) -> 5'(M). T/ien = e 2nir ^M (Wir) T {Pt0 rj P . 

Lemma 4.12 Fix ti > 1 1 < u 4 < ^ and < e < 1. Choose ?7i , T72 G <S(M) toii/t values in [0, 1], 

77l(i) = i IOrt(z[ ^K'2uK) md , ) = J M^G[ 2 L , 2uL j 

/1U \0 /or^[0,^) m j \0 /or^[0,f) 

and 1^1(01 <Ce-^l 1_e , (J 7 - 1 ^)! < Ce"^! 1- ' . Let rj P = r/i®n 2 . Thensuppnp C [0,^)x[0,f) = 
-Ra'.l ar *d the operator P G OPW has the following properties: 

a) The operator family 



^M uKk T^_ m _ y _L l PTuL l MK n _ uKk }^ ^ 



zs an 



l v ~ q -Riesz basis sequence for OPWw q {R 2 ) 



b) P G OPPF 1 ' 1 (i?ii- i L) and i/iere exisi functions di, d 2 : K — )■ Kq toii/i 

|P/0)| < 11/Hjifoo.oo di(x) and |P/(0l < d 2 (0, f G M 00 ' 00 ^), 

and di(x) < Ce-^l^l 1 "', d 2 (£) < Ce"^ 1- ' /or some C,7 > 0. 

Proof. The existence of 771, 772 satisfying the hypotheses stated above is established through mollifying 
characteristic functions. In fact, using constructions of Gevrey class functions, it has been shown 
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that for e, S > 0, there exists ip : M. — > Mq" and C, 7 > with suppc/? C [—6, 5], ftp = 1, < 
(7g-7l?l 1 6 I IihQ.'U IHor051 IDH98j . Note that the restriction to w subexponential in Theorem 



4.7 



is a 



consequence to the fact that there exist no compactly supported functions whose Fourier transforms 
decay exponentially (see references in |GP08j ). 
a) Due to Lemma |4.11 



being an Z~ ,l3 -Riesz basis for its closed linear span in m[^ ,uw (R 2 ) implies that 



r(l,l),(p,9) /'Tip 2 N 



1VI ^fn-uKk , , , 



is an i~ 9 -Riesz basis for its closed linear span in OPW2' 9 (M 2 ). 

6j As shown in the proof of Lemma 3.4 in [KP06J, we have \Pf(x)\ < |r/ 2 (— x)| II/IIm^^II^iIIm 1 ' 1 ; 
so we can choose d\{x) = \rf2{— %)\ WviWm 1 - 1 - 

Similarly, we can compute |P/(£)| < ||/||m°°>°° \\V2{£ ~ (") II Af 1 - 1 ■ W e claim that d.2(£) = \\t]2{€ ~ 
•)r/i(-)|| M i,i has the postulated subexponential decay. Recall that for g supported on [a, b], we have 



HsIImm < c II^IIl 1 where c depends only on the support size b — a (see Lemma 2.3 and [Oko09|). As 
t?2(£ — •)%(") * s compactly supported with uniform support size, we can compute 



d2(0 



■)m(-)\\M^<c\\T- 1 { m ^- •)%(•))! 



7/ 2 (£ - v)rji{v)e 



2irixv 



dv 



dx 



dx , 



(34) 



where 772(C) = 7/2 (~0- As the compact support of 771, r/ 2 together with IJ 7 //!^)! < Ce 7 '^' 1 " , \F 1 r/2(x)| < 
C e _ 7kl e imply that 771,772 are in the Gelfand-Shilov class S\zl [GS68J, we apply Proposition 3.12 in 
|GZ04] to conclude that VfT 2 7/i 6 <->ir|, that is, 

d 2 (0 < c J C^e-^'rt^ dx < Ce-^ |1_ \ 

□ 



Theorem 4.13 extends the main result in [Pfa08bJ to weighted and mixed Z p spaces with subexpo- 
nential weights. Both results generalize to infinite dimensions the fact that mx n matrices with m < n 
have a non-trivial kernel and, therefore, are not bounded below as operators acting on C n . 

Theorem 4.13 Let 1 < pi,P2,Qi,Q2 < 00, w\,W2 subexponential on l? d , that is, for some C, 7 > ; 
< P < 1, we have 



C-'e 



1 -7lMI 



< w\ (n) , W2 (n) < C e 



7IMI 



(35) 



If for M = (m fj ) : Cf 1 (Z 2d ) -> Cf 2 (Z 2d ), exists u > 1, K > 0, and 



with p<Ce^ H », /3 > /?, 



with 



\mfj\ < p(u\\j'\ 
then M has no bounded left inverses. 



3)5 «l|j'l 



> #0 
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Proof. Let v{n) = Ce^HI-. Note that l v ^ qx (l? d ) embeds continuously in Z^;°°(Z 2a! ) = l^ v {l? d ) 

and ^(Z 2 ^) = embeds continuously in l^ q2 (l? d ). Hence, it suffices to show that for all e > 

exists x G Z^(Z 2d ) with || x||;°° =1 and 1 1 Mx | |^i < e. For notational simplicity, we replace 2d by D in 
the following. 

First, observe that 

A Kl = e^" Yl K^e^ Y fc " 1 e~~^ -> as ^ -> oo. (36) 
Applying the integral criterion for sums, this would follow from 

R TOO /'OO ~ 

e7 Kf / x D-i e ^ y D ~ x dy dx ^ Q asKx -> oo. (37) 
For large substitution yields 

= if (D - 1)(1 -t } r(2^-D+2,7^) < if (Z? - 1)(1 -^(^) 2 ¥- D+1 e -V 
07 P /?7 

< hy D ~ 2 x 2D-l-P(D+l) e -7x/3 



where T denotes the upper incomplete Gamma function, and we used ya l\'l-y — > 1 as y — > 00 |AS92| . 



cbina vv c uocu ; 

The fact that (3 > (3 allows us to estimate for large x 

Hence, we can repeat the arguments above to the outer integral and the product with e 7 1 in (37) to 
obtain (37), and, hence, (36) holds. 

To continue our proof, we fix e > and note that (36) provides us with a K\ > K$ satisfying 
A Kl <{2 2D D 2 CC 2 e^ Nfi )- 1 e. 



As in |Pfa08bj . set N = A ^+ 1J and TV = \f ] + K x . Then ^ < N < A ^_+ 2 > implies 
XN > XKx + A + N and N > K x + f + 1 > K x + \f] = N. Hence, (2N + 1) D < (2N + 1) D so that 
the matrix M = (m.-/,-),, , n^» r : C^ 2Af+1 ^° — ?■ c( 2Ar+1 ) D has a nontrivial kernel. We now choose 
x G 1^ v (Ij D ) with ||x||/^ = 1, = if ||j||oo > N, and and Mx = where x is x restricted to the set 
{j : IUH00 < N} . 

By construction we have (Mx)ji = for ||j'||oo < N. To estimate (Mx)ji for ||j'||oo > AT> we 
fix K > K x and one of the 2£>(2(ff ] + if)) 5 " 1 indices j' G Z D with H/Hoo = [f ] + K. We have 
II Aj'Hoo > N + KX and XWfW^ - Wj]^ > KX > K for all j G iP with HjHoo < N. Therefore, using 
Holder's inequality for weighted /^-spaces, we obtain 

\( M x)j'\ = j Y m o'3 x \ ^ \\ x h? /v Yl v< j) \ m fj\ 

\\j\\oc<N \\j\\™<N 

< Ce^ P(A||J1» HIjIIdo) < C^"' Y P(\\j\U = 2 D DCe-< N "Y kD ~ 1 P( k )- 

llilU<Af Hillock k>K 
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Next, we compute 



|Mar||ji = J2vV')\(Mx) f \= £ v(j')\(Mx) f 



i'^ u lbl«,>Mrl+* 



< 2 D DCe^ £ E 

[|i'lloo>lTl+*i k >Wi'W°° 

< 2 D DCe* Nf> Yl 2D(2K) D - l C e^ K " £ k D ~ l p{k) 



K>\f]+K! k>K 



< 2 2D D 2 CC 2 e^ N " Yl K D - l e^ K "Y kD ' le ~^^ e - 

K >l^]+K! k>K 



□ 



Combining the results above, we can now proceed to prove Theorem 4.7 



Proof of Theorem J^.l As w is subexponential, there exists C, 7, e > with \w(x, £)| < C e^ x '^ c 



For this e > choose u, 771, 7/2, P, d±, and di as in Lemma |4.12| 

Define the synthesis operator E : C' 9 (Z 2 ) OPW% q {V M ) Q OPW$ q {M) as follows. For a = 
Wk,p} £ ^'(Z 2 ) write a^ p = ak : u+j for / € Z and < j < J and define 

_ J-i 

jB ( fJ ) = EE "W+i MuKkT^^PT^MK^^ (38) 
k,iez j=o 

with convergence in case p,q ^ 00 and weak-* convergence else. Since 

^ A^fTfc Tjl. si, PTul, Mk Kk > 

I A m a- 1 L n nftK J k,l,m,n& 

is an i~' 9 -Riesz basis for its closed linear span in OPWS' 9 (M 2 ), so is its subset 

I k^J+a 1 A" aPj ) k,leZ,0<j<J 

and -E 1 is bounded and bounded below . 



By Theorem 2.5, the Gabor system (g, a'Zxfe'Z) = {M^a'TifQ} is an ^~ 9 -frame for any a',b' > 
h a'b' < 1, and we conclude that the analysis map given by 

C g : Af™(R) -»■ ^ 9 (Z 2 ), /^((/.M^fe n)) (39) 

I AJ < J fc , 



is bounded and bounded below since u 2 K^j = u 4 k < 1. 

For simplicity of notation, set a = -fT and /3 = j^j. Fix / 6 M 00 ' 00 (]R) and consider the composition 

lP>i(Z 2 ) A OPW£' 9 (M) ^ M& 9 (R) ^ Z?f(Z 2 ) 

cr H> Per ^ Eag ^ { ( Eag, M u 2 ak ,T u 2p V g) } k , v . 

It is easily computed that the operator C^o^foE is represented — with respect to the canonical basis 
{<5(- — n)} n of Z~ 5 (Z 2 ) — by the bi-infinite matrix 

M = (m k >,i>,k,lJ+j) = (( MuakTht^^jPT-upu Mvj__ uak f , M u 2 ak ,T u 2p V g)J . 
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Setting 



di 



max Tk 
i=o,...,J-i - 



we observe 



-u/3j 



PT-ufnjMp 



fi-uak ■ 



\lTlk' ,1' ,k,U+j\ < (T u i3(lJ+j){T, 

< ll/l 

< ||/||m— - (U + j))) 



and 



(T uak M_ 



<Tu 2 /3i>9 ) 



-U/3U 



(PT_ uf 3u Mpj__ uak f )~, T u 2 ak ,M_ u ip V Q ) 



Uttfe 



/3J 



- u 2 ak' 



< (d 2 *fl)(ua(«A;'- &)). 

Observing that for appropriate C, 7, we have di *0(x)c?2 < C e^^H^'^^l 00 and 1 — 2e < 1 — e, 

□ 



allows us to apply Theorem 4.13 to Ai. This completes the proof. 



5 Outlook 

Recent results in operator identification with relevance to operator sampling include, for example, 
the identification of Multiple Input Multiple Output (MIMO) channels [Pfa08a] and an extension of 
operator sampling to irregular sampling sets [HP 09 . Some more fundamental questions concerning 
sampling and identification of operator Paley- Wiener spaces are still open. In the following we describe 
two such questions. 

Unbounded spreading domains with small Lebesgue measure. 



The extension of Theorem 4.6 to OPWw' q (M) with M unbounded but with Lebesgue measure less than 



one remains open. The following observations encourage tackling this question: 

1. Multiplication operators with not necessarily bandlimited symbol in L 2 (M 2 ) are clearly identifiable 
with identifier g = xr £ M°°'°°(M). Note that the characteristic function xr is the weak-* limit 
of T Y.nez S nT as T -> 0. Hence, the space OPW 00 ' 00 ({0} xl) is identifiable. 

2. Time-invariant operators with not necessarily compactly supported L 2 (M 2 ) impulse response are 
identifiable with identifier 5 which is the weak-* limit of Ylnez^nT as T — > 00. Consequently, the 
space OPW oo ' oo (]Rx{0}) is identifiable. 



3. In [KP06] it is shown that OPW(M) is identifiable if M is a possibly unbounded fundamental 
domain a lattice A in M 2 with A having density less than or equal to one. This result covers, 
for example, OPW({(t,u) : t> - 1, v<l, 2- ( - t+1 1<v<2- t }) as the unbounded set {{t,u) : 
t>-l, v<l, 2~(* +1 )<zv<2-*} is a fundamental domain of Z . 

The natural approach to construct identifiers for OPW(A) as weak-* limit of identifiers gjy for 



OPW(A n [-N, N] x [-N, N]) is difficult as the constants implied by x in (26) depend in a non-trivial 
matter on g^ = J2 n c n,N<5x n N , N G N, if the sequences {c nj jv} are not constant. 
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Generalizations to higher dimensions. 



As mentioned in Section |4j our proof of Theorem [46] hinges on the existence of identifiers in an analogous 
setup where the locally compact Abelian group E is replaced by an a ppropriate finite cyclic group of 
prime order Z p [KPR08, LPW05J. In fact, generalizing Theorem 4.6, to operators acting on L 2 (R d ) 



would be possible if the conclusions of Theorem 4.9 hold for sufficiently many composites n taking the 
place of prime p. In fact, in [KPR08 , we ask the following 



Question 5.1 Is it true that for all L £ N exists c G C L so that the vectors n(k, i)c, k,£ G 0, . . . , L — 1, 
defined by (ir(k,£)c)j = Cj-k e m T , k, £ = 0, . . . , L — 1, are in general linear position. 
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